Time resolved nanosecond vortex dynamics in high Tc superconducting films 
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We report on the observation of the rearrangement of vortices in the ns time scale following a fast 
magnetic field variation due to selective local heating of a high Tc superconducting film. Detailed 
simulations describing the measured photoinduced voltage signals along the film lead to a better 
understanding of vortex dynamics in fast regimes, where mean vortex velocities reach fO^ m/s. 



PACS numbers; 74.60. Ge, 74.60. Jg 

Vortex dynamics in type II superconductors continues 
to be a stimulating research area. Striking memory ef- 
fects and dynamical ordering of vortex matter in the low 
frequency domain have been observed recently |jl|], and 
may lead in the future to switching devices [^. However, 
the high frequency response of vortex matter is still an 
open question that needs to be addressed. 

In the past, the use of conventional techniques has 
limited the observations of magnetic properties to the 
low frequency domain. These techniques have long in- 
trinsic response and initial data-acquisition time&\, some- 
where between fractions of seconds and minutes Q . The 
study of short time vortex dynamics with novel non- 
conventional techniques is a requirement for the devel- 
opment of future fast devices. 

In this letter we report on the rearrangement of vor- 
tices in a high Tc superconducting sample over nanosec- 
ond time scales observed following a fast magnetic field 
variation due to partial heating. 

The sample is a GdBa2Cu307_a; film of thickness 
S = 300 nm on NdGaOa substrate. The sample was 
patterned (see Fig. |l|) using laser ablation in a strip 
(~ 10 X 2 mm^), and voltage contact pads were provided 
by gold sputtering, to reduce contact resistance below 
1 n. The Tc of the film is 91 K with a ATc = 1.5 K (10% 
- 90% criterion) determined by ac susceptibility. 

For the experimental observations we used a high 
bandwidth electronic detection method coupled with 
a pulsed magnetic field generator and a synchronized 
pulsed laser, based on a previously designed short-time 
magnetometer (~ 200 /iS ) Q]. The excimer laser pro- 
vided 45 ns FWHM pulses at 308 nm and uniformly il- 
luminated the target area. 

Briefly, the technique consists in the preparation of 
the thermomagnetic initial state by stabilizing the sam- 
ple temperature below the irreversibility line and apply- 
ing a transverse pulsed magnetic field Ha at a high rate 
{dH/dt ^ 1000 T/s) to reduce the uncertainty in the time 
origin. The laser pulse is triggered at a controlled delay, 
td, at different stages of the propagation of the magnetic 
field front. Flux penetration can be followed, as will be 
described below, by measuring the voltage signal across 
the pair of contacts provided in the film. The optical 
heating lowers the pinning force on the vortex structure 
(at the location where the laser spot was positioned) be- 



low that corresponding to the critical state at uniform 
temperature. This allows the fast penetration of vortices 
and triggers the flux redistribution in the whole sample. 
Flux motion is detected as a voltage pulse, superposed to 
signal induced by the application of the magnetic field. 

It has been shown that the photoinduced voltage is re- 
lated to the fast motion of vortices A thermal force 
cannot explain the results and instead a magnetic force 
drives flux motion. Having control of the time delay and 
the location of the illuminated area of the film, it is pos- 
sible to perform experiments that provide additional in- 
formation on fast vortex dynamics 
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FIG. 1. Schematics of experimental geometry and simula- 
tions. The mobile laser spot is shown in dashed lines. The 
measure loop ABCD has dimensions L and d. 

The model in study can be understood from Fig. 
Assuming homogeneity in the z direction and uniform 
applied transverse magnetic field in the x direction, the 
voltage reading is given by M: 



dt 
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where Ef, is the electric field at the contacts edge of the 
strip, Ba is the applied induction field, Bj is the field 
produced by the current distribution, L and d are the 
dimensions of the ABCD loop. We are considering B = 
fioH. The electric field Eb can be interpreted as given 
by the Josephson relation ^ = —v x B, where is an 
average vortex velocity and B the macroscopic magnetic 
induction. The first term in Eq. (|^) is a superposition 
of voltage pulses produced by individual vortices moving 
across the edge, each contributing in one flux quantum 
to the time integral of V{t). The last term in Eq. (|l|) is 
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the contribution from the change of the flux through the 
ABCD loop due to the return flux from the vortices. A 
moving vortex that enters the sample from the left will 
contribute only to this term. 

Assuming a phenomenological relationship between 
the electric field and the current density E = E{J, H,T), 
the basic features of the experimental results are well 
described by the quasi-stationary Maxwell equation [0 : 

dE _ fdHa S r dJ du \ 

'd^ ^^''"['^^2^ J_,dtV^j) ^ ^ 

Our experiment samples both the high and low cur- 
rent response of the material. In the high current regime 
above the critical state current density Jc, the resistivity 
takes the flux flow value p^j — pj^H/Hc2, where H is 
the local field, Hc2 is the upper critical one, and p„ is 
the normal state resistivity. In the low current regime, 
J < Jc, the resistivity is due to thermally activated creep 
and results in an electric field of the form E = Ec{J/Jc)^ 
pOf , where n >> 1. Observe that the resistivity at each 
point depends on the local magnetic field, which in turn 
is sensitive to the current distribution over the whole 
sample, leading to a highly nonlinear and nonlocal dy- 
namics. The temperature dependence of the critical cur- 
rent is given by MT) = J,o(l - T/T,)^/^ |ll|, where 
Ti = 86 K is the observed irreversibility temperature 
above which no photo-signal is detected. For the crit- 
ical field, in the temperature range under study we can 
approximate Hc2(T) = {T^-T) 1.8 x 10^ Am-^K'^ @. 

For the temperature independent parameters of these 
simulations we choose n — 20, SJco = 25 kA/m, pn/S = 
100 ri and E^ = 0.01 V/m. All measurements (and calcu- 
lations) were made at an initial temperature Tq = 73.2 K. 
The time-dependent temperature that results after local 
optical heating (see below) determines the time depen- 
dence of Jc and Hc2- 

To account for the experimental data we must assume 
that the crossover from one current regime to the other 
is smooth. We achieve this by interpolating both regimes 
with a cubic form in a small range (~ 5% Jc). 

Phenomenological constitutive relations of the kind we 
have discussed break down in the limit of very weak mag- 
netic fields, i.e., when we are dealing with regions of the 
sample where the magnetic flux has not yet penetrated. 
For this reason, we assume that there is a cut-off in the 
constitutive relation, such that E = whenever H < Hi, 
with Hi = 10~^Hc2- This regime is relevant to the early 
stages of the simulation. 

The resulting nonlinear nonlocal magnetic flux diffu- 
sion equation was integrated numerically by means of a 
time adaptive algorithm for the generally non-symmetric 
sheet current distribution 5J{y,t) [ |l3t . 

Figure || shows the voltage signal induced by the appli- 
cation of 7500 A/m. The rising edge was approximately 
exponential with a characteristic rise time t ~ 3.5 /is. 
Symbols indicate the measured voltage, and in full line 
we show the calculated voltage signal with an applied 



field modeled as an algebraic initial slope followed by an 
exponential rising edge. The lower signals, from bottom 
to top, correspond to the voltage on the edge of the sam- 
ple, and the difference with the full signal. Observe the 
delay in the onset of the edge electric field, associated 
with the formation of the characteristic current profile. 
Once this profile is formed, the model predicts the con- 
tacts voltage to be proportional to dHa/dt and indepen- 
dent of the sample parameters for the regime where the 
applied field changes faster than the induced one. 
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FIG. 2. Observed and simulated signal for the turning-on 
of the applied field. From bottom to top, the lines represent 
the contribution of the electric field at the edge of the film, the 
contribution due to the flux variation into the ABCD loop, 
and the total signal. We define t = at the maximum value 
of the measured signal. 
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FIG. 3. Voltage signals after optical heating of the sample. 
Open circles correspond to the illumination at the contacts 
side, and filled circles to illumination at the opposite one. 
Full lines show the simulated signals. The insets show the 
composition of these signals (see text). 

Illumination of the two thirds of the sample surface on 
the side of the contacts (right), 4 ps after the applica- 
tion of the field, produced the voltage pulse plotted in 
Fig. H (open circles); note the different scales. The local 
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temperature rise reduces the local critical current density 
below the flowing current density value (we neglect lat- 
eral thermal diffusion that would represent an increase 
in temperature in an additional small percentage of the 
non illuminated area during the duration of the pulse). 
This leads to a rapid reduction of sheet current density 
and the correlated penetration of magnetic flux. 

The curves plotted in the upper inset of the figure 
are the calculated voltage signals for this situation. The 
dashed curve corresponds to the electric field contribu- 
tion at the sample edge, the negative dotted curve is the 
term related to the magnetic flux variation through the 
ABCD loop, and the full line is the total calculated volt- 
age signal also shown in the main panel. 
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FIG. 4. Scanning of tlie film, varying tlie heated area. The 
circles show the time integrated signals. Open (filled) circles 
correspond to heating from the contacts (opposite to contacts) 
side. The full lines show the simulated results. The horizontal 
dashed line shows the value obtained without heating. The 
vertical line shows the heated area chosen for other figures. 

The filled circles show the corresponding result when a 
similar sample surface is illuminated at the edge opposite 
to the contacts (the left one). Notice that the measured 
voltage changes sign. This result is well described by the 
calculations (lower inset), where the flux time derivative 
is almost the only contribution to the total voltage . In- 
deed, almost no vortices move across the contacts edge 
as they are magnetically forced to enter the fllm from the 
opposite direction. 

To test the model further, the laser spot was scanned 
across the film, repeating the experiment at a fixed tem- 
perature and at a fixed delay time td- The measured 
voltage signals were time integrated {Vint) and are plot- 
ted in Fig. ^ (circles) as a function of the percentage of 
heated area of the strip. While the upper curve corre- 
sponds to heating the contacts side of the sample, the 
lower curve corresponds to heating the opposite side. 

The temperature time dependence at the laser spot for 
these simulations is: 
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where Tq and Ti are the initial and the irreversibility 
temperatures, and t is the time elapsed since the laser 
pulse is triggered. The parameter values are N — 4, C = 
15 ns, M = 4.3 and A = 2.8. The essential feature of this 
curve is the asymmetry between the fast heating and the 
slow cooling of the illuminated area. This behavior has 
been observed in simulations of heat diffusion for realistic 
values of diffusion coefficients of sample and substrate, 
and sample thicknes s |14[] . The particular values of the 
coefficients in Eq. (Q) follow from a global fit to the 
experimental data. Observe that in our simulations the 
temperature rises from the initial value Tq — 73.2 K to 
Tmax = 90 K, which is higher than the irreversibility 
temperature T!^ = 86 K. 

From the numerical solution of the evolution equation 
(Eq. (||)) with the constitutive relation already discussed 
we have reconstructed the time evolution of the sheet cur- 
rent and the magnetic and electric fields inside the sam- 
ple. Also, assuming the Josephson relation E = —v x B, 
we have obtained the mean velocity field of the vortices 
within the sample. 
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FIG. 5. Space distribution of current density (upper), mag- 
netic field (middle) and vortices velocities (lower) for different 
times since the start of the laser pulse. The inset in the lower 
panel details the symmetry of the initial condition. 

In Fig. Ijwe show the calculated sheet current distribu- 
tion (upper panel) at different times since the triggering 
of the laser pulse. The vertical line shows the border be- 
tween cold (left) and heated (right) areas. We also show 
the magnetic field (middle panel) and the velocity field 
(lower panel) at the corresponding times. Observe that 
a positive sign in this last plot represents left moving 
vortices, while the applied field points in the negative x 
direction. Note that velocities after heating reach values 
of order 10** m/s and before heating 10^ m/s. To appre- 
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ciate the different scales we have inserted in the lower 
figure the plot of the initial velocities with a higher res- 
olution. At td the velocities are odd because the applied 
field is homogeneous. The calculated mean velocities are 
consistent with earlier observations of propagation of flux 
fronts 

The rearrangement of the vortices is triggered by the 
fast decay in the critical current produced by heating. 
In fact, the critical sheet current drops to zero, since 
the temperature goes over the irreversibility line. While 
magnetic field penetration is easier from the right edge of 
the strip, some vortices enter the illuminated area from 
the left, meaning that we must have vortex motion in the 
cold area as well. 

Observe that a possible temperature dependence in the 
creep exponent n would not affect our results, since in the 
cold zone temperature is essentially constant, and in the 
hot zone the sheet current is over the critical one and so 
in the flux flow regime. 

Our simulations also display two slower phenomena, 
not shown in Fig. ^. These are the relaxation towards 
an homogeneous distribution of vortices in the hot zone, 
with a characteristic time of the order of the /iS, and then 
the change in magnetization in response to the changing 
applied field, which is even slower. 

In conclusion, we have developed a non-conventional 
technique whereby an abrupt change in magnetic field 
on the sample is accomplished not by varying the ap- 
plied field, but by heating the sample with a laser pulse. 
By scanning the laser pulse across the sample, delaying 
the application of the pulse with respect to the applica- 
tion of the external magnetic field, and by heating the 
sample both on and across the contacts' side, we have 
been able to observe different aspects of vortex dynamics 
on nanosecond time scales, and to discriminate between 
the different contributions to the measured signal. By 
means of numerical simulations, we have been able to re- 
construct the evolution of the current and vortex motion 
within the sample in these fast time scales, thus opening 
a new time window for the study of vortex dynamics in 
high Tc superconductors. 
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